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1. The scaled ABS class. 

ABS methods have been introduced by Abaffy, Broyden and Spedicato (1984), in a paper 
where the solution of linear algebraic equations was considered via the so called ABS un- 
sealed or basic class. That class was later generalized to the so called scaled ABS class and 
then applied also to the solution of linear least squares, nonlinear algebraic equations and 
optimization problems. For a full presentation of the main results in the ABS field see the 
monographs of Abaffy and Spedicato (1989) and of Zhang, Xia and Feng (1999). There are 
presently almost 400 papers dealing with ABS methods, see for a partial list the bibliography 
in the Appendix 2. In this review we restrict our attention to ABS methods for linear de- 
termined or underdetermined algebraic systems and their application to some optimization 
problems. We begin by giving our notation and the steps of the scaled ABS class. 

Consider the following general linear system (determined or underdetermined), where the 
rank of A G Jl m > n [ s arbitrary 



Ax = b x G R n , be R m , m<n 



or 



cij x — hi = 0, i = 1, 



, m 



where 



A 



T 



(1) 
(2) 

(3) 



Then the scaled ABS class is defined by the following procedure: 

1. Let x\ G R n be arbitrary, let H\ G R n ' n be nonsingular arbitrary, let v\ G R m be 
arbitrary nonzero, set i = 1. 
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2. Compute the residual vector = Axi — b. If = stop (xj solves the system), 
otherwise compute Sj = HiA T Vi. If Sj 7^ 0, go to (3). If Sj = and r = vJYj = 0, set 
Xi + \ = Xi, H i+ i = ilj and go to (6). Otherwise stop, the system has no solution. 

3. Compute the search vector p, t by 

Pi = Hfz t (4) 

with Zi G R n arbitrary save that 

vjAHfzi + (5) 



4. Update the estimate of the solution by 

Xi+\ = Xi- aipi (6) 

where the stepsize on is given by 

a-i = vfn/vfApi. (7) 

5. Update the matrix Hi by 

H i+1 = Hi- H l A T v i wfH l /wfH l A T v l (8) 
where Wi G i? n is arbitrary save that 

w?HiA T Vl + 0. (9) 



6. If i = m, then stop (x m+ i solves the system), otherwise give Vi+i G R m arbitrary 
linearly independent from v±, . . . ,Vi, increase i by one and go to (2). 

Matrices Hi are generalizations of oblique projectors. They were first introduced in literature 
apparently by Wedderburn (1934) and have been used in several papers of Egervary that in 
fact introduced ABS methods ante litteram. They have been denoted as Abaffian matrices at 
the First International Conference on ABS Methods (Luoyang, China, 1991) and this name 
will be used here. 

The above recursions define a class of algorithms. A particular method is obtained by a 
choice of the parameters H 1 , Vi, ,Zi, Wi. The basic A BS class is obtained by setting vi = e^, 
being the z-th unit vector in R m . Parameters Wi,Zi,Hi were introduced respectively by 
Abaffy, Broyden and Spedicato, from whose initials the class name is obtained. One can 
show that the scaled ABS class is a complete realization of the so called Petrov-Galerkin 
principle for solving a linear system, where the iteration has the form Xj+i = x% — onpi with 
Qj, pi such that r^Vj = 0, j = 1, . . . , i, where vectors Vj are arbitrary linearly independent. 
It appears that all methods proposed in the literature for solving a linear system in a finite 
number of steps are special cases of the ABS class. Indeed even the Quasi-Newton methods 
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of the Broyden's class for nonlinear systems, which are known to have finite termination 
under mild conditions in no more that 2n steps on linear systems, have been shown recently 
to be imbeddable in the ABS class, see Adib, Mahdavi Amiri and Spedicato (1999), in the 
sense that the approximation to the solution of index 2i — 1 generated by Broyden method is 
identical with the approximation of index i of a certain ABS algorithm, when both algorithms 
start from a same initial point. 

The recursions of the scaled ABS class can be derived in different ways. For instance, let 
us consider the basic ABS class, which has the property that the approximation £1+1 solves 
the first i equations of the system Ax = b. This property is obtained by a suitable choice of 
the stepsize and the search vector in the update Xj+i = Xi — ctiPi. Choosing aj by formula 
(7) makes the i— th equation a[x = bi to be satisfied by Xj+i. Choosing p, L by (4) with Hi 
updated by (8) makes the previous i — 1 equations to remain satisfied in the new point, since 
a simple induction argument shows that this is true if HiCij = for j < i, a property that 
characterizes update (8). Conditions (5) and (9) guarantee that the process is well-defined 
(no division by zero). Condition (5) has a precise geometrical justification, implying that the 
i— th search vector intersects the linear variety containing all solutions of the i— th equation 
at a finite point (the vector being orthogonal to such a variety). The recursions of the 
scaled ABS class are obtained from those of the basic ABS class by applying them on the 
system V T Ax = V T b, where V = (vi,..., v m ), which is equivalent in the sense of possessing 
the same set of solutions, if V is nonsingular. 

Referring for proofs to the monograph of Abaffy and Spedicato (1989), we recall now 
some properties of the scaled ABS class. For simplicity we assume that rank(A) = m. 

• Define Vi = (v\, . . . ,Vi), Wi = (wi, . . . ,Wi). Then Hi + \A T Vi = 0, Hf +l Wi = 0, hence 
vectors A T Vj,Wj, j = 1, . . . ,i span respectively Null(Hi + i) and Null(Hi + i) T . 

• Vectors HiA T Vi, Hfwi are zero iff Oi, Wi respectively linearly depend on a±, . . . , Oj_i, 

Wi, . . . ,Wi-!. 

• Define Pj = (p±, . . . ,pi) and Ai = (a±, . . . , ctj). Then we have the implicit factorization 
Vi Aj Pi = Li, with Li nonsingular lower triangular. Hence for m = n we get a 
semiexplicit factorization of the inverse of A, i.e. with P = P n , V = V n , L = L n 

A' 1 = PL~ l V T . (10) 

For several choices of V, matrix L is diagonal, hence formula (10) gives an explicit 
factorization of the inverse, a bonus freely provided by the ABS procedure. One can 
moreover show that all possible factorizations of a matrix can be obtained by suitable 
choices of the ABS parameters, another completeness result. 

• Define Si and Pj as Si = (s\, . . . , Sj), Pj = (ri, . . . , j-j), where Sj = HiA T Vi, rj = Hjuti. 
Then P«+i = H\ — SiRj and vectors Sj, r,- L and p, L can be built by a generalized Gram- 
Schmidt iteration using 2i vectors. This alternative formulation is computationally 
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convenient in term of storage when m < n/2. There is also an alternative formulation 
using n — i vectors. 

• One can write the Abafnan as follows, in term of the parameter matrices 



Setting V = V m , W = W m , one can prove that Hi, V, W are admissible (i.e. con- 
dition (9) is satisfied) iff matrix Q = V T AHfW is strongly nonsingular (i.e. Q is LU 
factorizable) . 

• If m < n the general solution of Ax = b can be written as follow, with q € R n arbitrary 



It was observed that the scaled ABS class can be obtained by applying the unsealed 
ABS class to the scaled system V T Ax = V T b. This scaled system can be viewed also as 
a left preconditioned system. Hence the scaled ABS class includes also all possible left 
preconditioners, which are defined in the ABS context implicitly and dynamically. 

2. Some subclasses and methods 

Here we consider some subclasses and special methods in the scaled ABS class. 

(a) The conjugate direction subclass. This subclass is given by the choice Vi = pi and is well 

defined under the sufficient but not necessary condition that A is symmetric positive 
definite. All methods in this subclass have search directions that are ^4-conjugate 
and conversely any method that generates A-conjugate directions is a member of the 
subclass. If xi = then Xj+i is the vector of least energy norm (i.e. ^4-weighted 
Euclidean norm) in the linear variety xi + Span(pi, . . . ,pA and moreover the solution is 
approached monotonically in such a norm. The subclass contains the ABS reformulation 
of the Choleski algorithm and of the methods of Hestenes-Stiefel and of Lanczos. 

(b) The orthogonally scaled subclass. This subclass is obtained by the choice V{ = Api 

and is well defined if A has full column rank; it is also well defined even if m > 
n, in which case the iteration can be performed for n steps and will stop in a least 
squares solution. The subclass generates scaling vectors that are orthogonal and search 
directions that are yl T yl-conjugate. If xi = then xi + i minimizes the Euclidean norm 
of the residual in the linear variety xi + Span{jpi, . . . ,p{), implying that the solution is 
monotonically approached from below in the residual norm. This subclass contains the 
ABS formulation of the QR and GMRES methods. 

(c) The optimally stable subclass. This subclass is obtained by the choice vi = A~ T pi, the 

inverse disappearing in the actual recursion formulas. The class is well defined without 



H i+ i = Hi- HiA^iiWfH^V^WfHi. 



(11) 



x = x m+ i + H^ +1 q 



(12) 
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conditions on the matrix (assumed to be square nonsingular, a condition that can be 
relaxed). The search vectors in the subclass are orthogonal. If x\ = then Xi+± is 
the vector of least Euclidean norm that solves the first i equations, implying that the 
solution is monotonically approached from below in the Euclidean norm. The name 
of the subclass derives from the fact that a condition is satisfied which guarantees the 
least propagation of a single error introduced in the sequence X{. This subclass contains 
the Huang algorithm, see the following paragraph, and the ABS version of the Craig 
algorithm. 

(d) The Huang method. This method belongs to both the basic ABS class and the optimally 

scaled subclass. It can be obtained by the parameter choices H\ = I , Z{ = Wi = cij, V{ = 
ej. A mathematically equivalent, in the sense of generating the same iterates in exact 
arithmetic, but numerically more stable form is the modified Huang method, based upon 
formulas pi = Hi(Hiai) and = Hi — Pipf /pfpi- The Huang method generates 

orthogonal search vectors. If x\ = then Xi+i is the solution of least Euclidean norm 
of the first i equations and the solution x + of least Euclidean norm of the system is 
approached monotonically and from below from sequence Xj. Zhang (1995a) has shown 
that using the active set strategy of Goldfarb and Idnani (1983) the Huang algorithm 
can solve a system of general linear inequalities where in a finite number of steps it 
either establishes that no feasible point exists or it finds the solution of least Euclidean 
norm. If one takes as initial matrix a symmetric positive definite matrix then all 
above properties remain valid with respect to the -B-weighted Euclidean norm. 

Extensive numerical experiments, see Spedicato, Bodon and Luksan (2000) have shown 
that a certain version of the modified Huang algorithm is very accurate, being compa- 
rable on extremely ill conditioned problems with the LAPACK implementations of the 
rank revealing QR factorization or the SVD factorization. It provides usually the same 
estimate of the numerical rank as these LAPACK codes. Its computational times are 
up to over one hundred times lower on problems with 2000 variables. 

(e) The implicit LU method. This method is obtained by the choices H\ = I, Zi = Wi = 

Vi = ei, which are well defined iff A is strongly nonsingular (otherwise column pivoting is 
needed, or row pivoting if m = n). The parameter choices imply the following structure 
for the Abaffian, with K { € R 71 '^ 



H^= ^ Q • (13) 
K T 

1 n—t 

From (13) we get that the last n — i components of pi are zero, while the i-th component 
equals one, hence, if x\ = 0, the sequence Xi consists of special basic vectors. Moreover 
Pi is upper triangular with units on diagonal, hence the associated implicit factorization 
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is of the LU type, justifying the name. If m = n, then n 3 /3 multiplications are required, 
plus lower order terms. Main storage is for Ki, which is at most n 2 /4. Hence the 
algorithm can be considered more efficient than classical Gaussian elimination or LU 
factorization, since it requires the same number of multiplications but has a lower 
storage need (in the case that A cannot be overwritten). 

(f) The implicit LX method is given by choices H\ = I,Vi = ej, Z{ = wi = e^, with ki, 
1 < ki < n such that 

elHiCH + 0. (14) 

If A has full rank, then there exists at least one index fej such that (14) is satisfied. For 
stability we usually choose ki to maximize | HiOi |. 

Let N = (1, 2, . . . , n), let B i = (k 1 ,..., h) and let N t = N\Bi. Then we have: 

• ki e Ni-i 

• Hf +1 e k = for k € B { 

• Vector pi has n — i zero components; the fcj-th component equals one 

• If X\ = 0, then Xi + \ is a basic type solution of first % equations 

• Columns of ilj+i of index k G Ni are the unit vectors e^; columns of index k G Bi 
have zero components in j'-th position, with j £ Bi. 

• At i-th step i(n — i) multiplications are required to compute H-iai and i{n — i) to 
update the non trivial part of Hi. Hence the operation cost is the same as for the 
implicit LU algorithm (i.e. n 3 /3), but pivoting is not required. 

• Storage is again the same as for the implicit LU algorithm, i.e. at most n 2 /4. 

It follows from the above that the implicit LX algorithm can be considered as the most 
efficient general solver for a linear system not of the Strassen type, since it has the same 
operation cost of Gaussian elimination, requires less storage (when A is not overwritten) 
and requires no condition on the matrix. It should also be remarked that numerical 
experiments by Mirnia (1996) have shown that it provides a very accurate estimate 
of the solution, usually better than by the implicit LU methods or the classical LU 
factorization as implemented by MATLAB or LAPACK. 

3. Other ABS algorithms and numerical experiments 

For most algorithms proposed in the literature their ABS reformulation has been obtained, 
leading, in view of the several existing alternative formulations of the ABS linear algebra, to 
new formulations in many cases with better computation cost, storage requirement, numerical 
stability and possibility of vectorization and parallelization. The reprojection technique used 
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in the modified Huang algorithm with success, based on the identities R; L q = H^H^), Hfq = 
Hf(Hfq), valid for any vector q if Hi = I, has been used to improve substantially the stability 
of the methods of Hestenes-Stiefel and of Craig. Moreover the classical iterative refinement 
can be reformulated in the ABS context using the previous search vectors. Experiments 
by Mirnia (1996) have shown that such refinement coupled with residual computation in 
double precision provides usually a fully accurate single precision solution in just one or two 
refinement steps. 

For structured systems one can develop special ABS methods exploiting the problem 
structure; often the Abaffian matrix reflects the structure of A. For instance if A is banded, 
so is the Abaffian generated by implicit LU, implicit QR and Huang methods, albeit with a 
larger band width. If A is symmetric positive definite with a Nested Dissection structure, so 
is the Abaffian generated by the implcit Choleski algorithm, whose storage turns out to be, 
for n going to infinity, infinitesimal with respect to the storage required for the corresponding 
special implementation of the classical Choleski method. For matrices having the KT (Kuhn- 
Tucker) structure, large classes of ABS methods have been derived, many of them having up 
to one order less storage or number of operations than the classical methods based upon LU 
or QR factorization or the Aasen process, see Spedicato, Chen and Bodon (1996). Moreover 
testing of some of such ABS methods has also shown that they are usually more accurate, 
see Spedicato, Bodon and Luksan (2000). 

For matrices with a general sparsity pattern little is known about the fill-in structure of 
the Abaffian. We have however found that by clever use of BLAS4 routines the performance 
of the full version is greatly improved and becomes competitive with the performance of codes 
like MA28 minimizing the fill-in in LU factors, at least for n not too large. 

There are also block formulations of the ABS procedure, which do not lead to accuracy 
deterioration and provide reduction in number of operations and, in particular, faster im- 
plementations on vector or parallel computers. The special case of the two-block has been 
investigated by Adib and Mahdavi-Amiri (2000), with good numerical results. 

Finally ABS methods can be used to generate large classes of iterative methods for linear 
systems via two approaches. The first is based upon restarting the iteration after k < m 
steps and using the vector formulation instead of the full Abaffians, so that the storage is 
only of order 2k vectors. The second approach is based upon the truncation idea. The 
vector formulation contains summations, at the i— th step from 1 to i. Thus we only keep p 
terms in such summations, typically the last p terms, so that storage is of order p vectors. 
Global convergence at a Q-linear rate has been established for large subclasses of the scaled 
ABS class, including the conjugate direction, the orthogonally scaled and the optimally stable 
subclasses, under mild conditions on the parameter choices, following the technique of Dennis 
and Turner (1987). The obtained class of iterative methods includes such classical methods 
as the Krylov space methods, the Gauss-Seidel, the De La Garza and Kackmartz methods, 
see Spedicato and Li (2000). 

A project named ABSPACK is now under development aiming at producing documented 
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software for solving linear and nonlinear algebraic equations and some optimization problems 
using ABS methods. In the course of 1999 algorithms have been implemented for solving full 
linear underdetermined, determined and overdetermined (in the least squares sense) and KT 
systems. Comparisons have been made with other commercial packages, including LINPACK, 
LAPACK, MATLAB, UFO, showing that the produced codes are competitive in terms of 
both accuracy and timing (better results in timing are expected with the use of the block 
formulations, currently under development). In the Appendix some results are given, taken 
from Spedicato, Bodon and Luksan (2000). 

4. ABS solution of linear Diophantine equations 

The ABS algorithm has been recently applied to solve linear Diophantine systems, where it 
determines if the system has an integer solution, computes a particular solution and provides 
a representation of all integer solutions. The ABS algorithm is a generalization of a method 
proposed by Egervary (1955) for the particular case of a homogeneous system. 

Let Z be the set of all integers and consider the Diophantine linear system of equations 

Ax = b, x G Z n , A G Z mxn , beZ m , m<n. (15) 

We recall some results from number theory. Let a and b be integers. If there is integer 
7 so that b = 7a, then we say that a divides b and write a\b, otherwise we write a/f b. If 
a±,... ,a n are integers, not all being zero, then the greatest common divisor (gcd) of these 
numbers is the greatest positive integer 5 which divides all a^, i = 1, . . . , n and we write 
5 = gcd{a\ , . . . , a n ) . We note that 5 > 1 and that 5 can be written as an integer linear 
combination of the Oj, i.e. 5 = z T a for some z G R n . One can show that 5 is the least positive 
integer for which the equation a\X\ + . . . + a n x n = 5 has an integer solution. Now 8 plays a 
main role in the following 

Fundamental Theorem of the Linear Diophantine Equation 

Let a±, . . . , a n and b be integer numbers. Then the Diophantine linear equation a\X\ + . . . + 
a n x n = b has integer solutions if and only if gcd(a±, . . . ,a n )\ b. In such a case if n > 1 then 
there is an infinite number of integer solutions. 

In order to find the general integer solution of the Diophantine equation a\X\ + . . . + a n x n = 6, 
the main step is to solve a±x± + . . . + a n x n = 5, where 5 = gcd(a\ , . . . , a n ) , for a special integer 
solution. There exist several algorithms for this problem. The basic step is the computation 
of 5 and z, often done using the algorithm of Rosser (1941), which avoids a too rapid growth 
of the intermediate integers, and which terminates in polynomial time, as shown by Schrijver 
(1986). The scaled ABS algorithm can be applied to Diophantine equations via a special 
choice of its parameters, originating from the following considerations and Theorems. 

Suppose Xi is an integer vector (we begin with x\ an integer vector). Since Xj+i = Xi—ctiPi, 
then Xi + \ is integer if «j and pi are integers. If vfApi\(vfri) then on is an integer. If Hi and 
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Zj are respectively an integer matrix and an integer vector, then pi = Hj Zi is also an integer 
vector. Assume Hi is an integer matrix. From (6), if vf AHfwi divides all the components 
of HiA T ViwJ Hi, then H i+ i is an integer matrix. 

Conditions for the existence and determination of the general solution of the Diophan- 
tine system are given in the following theorems, generalizing the Fundamental Theorem, see 
Esmaeili, Mahdavi-Amiri and Spedicato (1999). 

Theorem 1 

Let A be full rank and suppose that the Diophantine system (15) is integerly solvable. Con- 
sider the Abaffians generated by the scaled ABS algorithm with the parameter choices: H\ is 
unimodular (i.e. both H\ and H^ 1 are integer matrices); for i = 1, . . . ,m, Wi is such that 
wj HiA T Vi = Si, Si = gcd(HiA T Vi) . Then the following properties are true: 

(a) the Abaffians generated by the algorithm are well-defined and are integer matrices 

(b) if y is a special integer solution of the first i equations, then any integer solution x of 

such equations can be written as x = y + Hf +1 q for some integer vector q. 

Theorem 2 

Let A be full rank and consider the sequence of matrices Hi generated by the scaled ABS 
algorithm with parameter choices as in Theorem 1. Let x\ in the scaled ABS algorithm be an 
arbitrary integer vector and let z,- L be such that zj HiA T V{ = gcd(HiA T Vi). Then system (15) 
has integer solutions iff gcd(HiA T Vi) divides vjri for i = 1, . . . , m. 

The scaled ABS algorithm for Diophantine equations is now stated. 

The ABS Algorithm for Diophantine Linear Equations 

(1) Choose x\ € Z n , arbitrary, H\ <G Z nxn , arbitrary unimodular. Let % = 1. 

(2) Compute r, = vjri and Sj = HiA T Vi. 

(3) If (si = and Tj = 0) then let xi + \ = Xi, Hi + \ = Hi, r,i + \ = r,- L and go to step (5) 
(the ith equation is redundant). If (si = and Tj 7^ 0) then Stop (the iih equation and 
hence the system is incompatible) . 

(4) {s{ 7^ 0} Compute Si = gcd(si) and pi = Hf Zi, where Zi G Z n is an arbitrary integer 
vector satisfying zj Si = Si. If Si /frj then Stop (the system is integerly inconsistent), 

H A T v w T H 

else Compute oti = n/Si, let x i+ i = Xi-onpi and update Hi by H i+1 = Hi- \ < „\t„ 1 
where Wi € R n is an arbitrary integer vector satisfying wfsi = Si. 

(5) If i = m then Stop (x m +i is a solution) else let i = i + 1 and go to step (2). 

It follows from Theorem 2 that if there exists a solution for system (15), then x = x m+ i + 
Hm+iQi with arbitrary q € Z n , provides all solutions of (15). 
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5. Application of ABS methods to unconstrained optimization 

A direct application of ABS methods to the problem 

minf = f(x), x G R n 

can be made via the iteration Xj+i = Xi — diHf Zi, where Hi is reset after n or less steps 
and Zi is chosen to satisfy the descent condition, i.e. gj Hj Zi > 0, with gi the gradient of 
j{x) in xi. If fix) is quadratic with Hessian A, we may identify matrix A in update (8) with 
such Hessian. Defining x' = x% — (3vi with (3 a small scalar and noticing that on quadratic 
functions one has g' = g — (3Avi, update (8) can be written without use of A as 

H i+1 = Hi- HiyiwJ Hi/wf Hiyi (16) 
with Hi = g' — gi, while the stepsize takes the form 

ai = vjn/pjyi (17). 

The above defined class has termination in at most n steps on quadratic functions and on 
general functions has local convergence at n-step Q-quadratic rate. It is related to a class of so 
called projection methods defined by Pschenichny and Danilin (1978). Numerical experience 
with such methods is not yet available. 



6. Derivation of general Quasi-Newton updates 

ABS methods allow to determine the general solution of the Quasi-Newton equation, even 
in presence of additional linear conditions as symmetry and sparsity. Moreover they provide 
insight into the problem of keeping positive definiteness. Notice that while the general solution 
of the Quasi-Newton equation was known since a paper of Adachi (1971), the ABS explicit 
solution for the problem with sparsity and symmetry was not obtained before. 

Let us consider the Quasi-Newton equation defining a new approximation to the Jacobian 
(for a nonlinear system, g(x) = 0) or to the Hessian (for an optimization problem, min f(x), 
with g(x) the gradient of f(x)), here given in the transpose form 

d T B' = y T (18) 

where d = x' — x,y = g' — g. We can notice that (18) is a set of n underdetermined linear 
systems, each one having just one single equation and differing only in the right hand side. 
Thus one step of the ABS procedure is enough to get a special solution and a representation 
of all solutions via formula (12). By suitable choice of the parameters the general solution 
can be put in the following form, see Spedicato and Xia (1993) 

B' = B — s(B T d - y) T /d T s + (I - sd T /d T s)Q (19) 
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with Q 6 R n ' n arbitrary and s € R n arbitrary save that s T d / 0. Formula (19) is the same 
formula derived by Adachi (1971). 

The condition that some elements of B' are zero or have a constant value and the condition 
of symmetry can be written as the following additional linear equations, where b[ is the i-th 
column of B' 

( b 'i) T ej = Vij 

where rjij = implies sparsity, rjij = costant implies constant value and rjij = rjji implies 
symmetry. The above conditions can be dealt with by the ABS algorithm, see Spedicato and 
Zhao (1993), obtaining an explicit solution, by successive columns in the case of symmetry. 
By requiring that the first n — 1 diagonal elements are sufficiently large, then we obtain 
formulas where B' is almost positive definite or almost diagonally dominant (i.e. its (n — 1)- 
th principal submatrix is positive definite or diagonally dominant). It is not possible in 
general to have B' positive definite, since the last column is completely defined by the n — 1 
symmetry conditions and the Quasi-Newton equation. 

7. The simplex method via the implicit LX algorithm 

The implicit LX algorithm has a natural application in reformulating the simplex method for 
the LP problem in standard form: 

T 

mm c x 

subject to 

Ax = b, x > 0. 

The suitability of the implicit LX algorithm derives from the fact that the algorithm, when 
started with x\ = 0, generates basic type vectors Xj+i, which are vertices of the polytope 
defined by the constraints of the LP problem if the components not identically zero are 
nonnegative. 

The basic step of the simplex method is moving from one vertex to another according 
to certain rules and reducing at each step the value of e x. The direction of movement is 
obtained by solving a linear system, whose coefficient matrix Ab, the basic matrix, is defined 
by m linearly independent columns of A, called the basic columns. This system is usually 
solved via the LU factorization method or also via the QR factorization method. The new 
vertex is associated with a new basic matrix A' B , obtained by substituting one of the columns 
of Ab with a column of the matrix An that comprises the columns of A not belonging to 
Ab- One has then to solve a new system where just one column has been changed. If the LU 
factorization method is used, then the most efficient way to recompute the modified factors 
is the steepest edge method of Forrest and Goldfarb (1992), that requires m? multiplications. 
This implementation of the simplex method has the following storage requirement: m 2 for 
the LU factors and mn for the matrix A, that has to be kept to provide the columns needed 
for the exchanges. 
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The reformulation of the simplex method via the implicit LX method, due to Xia (1995c), 
Zhang and Xia (1995), Spedicato, Xia and Zhang (1997), Feng et al. (1997), exploits the 
fact that in the implicit LX method the exchange of the j-th column in A B with the k- 
th column in An corresponds to exchanging previously chosen parameters Zj = Wj = ej B 
with new parameters zt = Wk = e^ B . In term of the implied modification of the Abaffian 
this operation is a special case of a general rank-one modification of the parameter matrix 
W = (wk-L, ■ ■ ■ , Wkm)- The modified Abaffian can be efficiently evaluated using a general 
formula of Zhang (1995d), without explicit use of the fc-th column of An. Moreover all 
information needed to build the search direction (the polytope edge) and to implement the 
(implicit) column exchange is contained in the Abaffian matrix. Thus there is no need to 
keep the matrix A. Hence storage requirement is only that needed in the construction and 
keeping of the matrix H m+ ±, i.e. respectively n 2 /4 and n(n — m). For m close to n storage, it 
is about 8 times less than in the Forrest and Goldfarb implementation of the LU method. For 
small m storage is higher, but there is an alternative formulation of the implicit LX method 
that has a similar storage, see Spedicato and Xia (1999). 

We now give the main formulas for the simplex method in the classical and in the ABS 
formulation. The column in An that substitutes a column in As is usually taken as the 
column with least relative cost. In the ABS approach this corresponds to minimize with 
respect to i G N m the scalar r\i = c T H T ei. Let N* be the index so chosen. The column in Ab 
to be exchanged is often chosen with the criterion of the least edge displacement that keeps 
nonnegative the basic variables. Defining uii = x T ei/ ej H T e^* with x the current vertex, 
then the above criterion is equivalent to minimize uj{ with respect to the set of indices i G B m 
such that 

efH T e N * > (20) 

Notice that H t cn* / and that an index i such that (20) be satisfied always exists, unless 
£ is a solution. 

The update of the Abaffian after the exchange of the obtained unit vectors is given by 
the following simple rank-one correction 

H' = H- (He B * - e B ')e£.H/e£.HeB* (21) 

The displacement vector d, classically obtained by solving system A B d = — Acn*, is obtained 
at no cost by d = H 7 +1 eN* • The relative cost vector, classically given by r = c — A 7 A^cb, 
with cb comprising the components of c with indices corresponding to those of the basic 
columns, is given by formula r = H m+ \c. 

It is easily seen that update (21) requires no more than m{n — m) multiplications. Cost 
is highest for m = n/2 and gets smaller as m gets smaller or closer to n. In the method 
of Forrest and Goldfarb m? multiplications are needed to update the LU factors and again 
m? to solve the system. Hence the ABS approach via formula (21) is faster for m > n/3. 
For m < n/3 similar costs are obtained using the alternative formulation of the implicit LX 
algorithm described in Spedicato and Xia (1999). For m very close to n the advantage of the 
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ABS formulation is essentially of one order and such that no need is seen to develop formulas 
for the sparse case. 

There is an ABS generalization of the simplex method based upon a modification of the 
Huang algorithm which is started by a certain singular matrix, see Zhang (1997). In such 
generalization the solution is approached via a sequence of points lying of the faces of the 
polytope. If one of such points happens to be a vertex, then all successive iterates are vertices 
and the standard simplex method is reobtained. 

8. ABS unification of feasible direction methods for linearly constrained mini- 
mization 

ABS methods allow a unification of feasible direction methods for linearly constrained mini- 
mization, of which the LP problem is a special case. Let us first consider the problem with 
only equality constraints 

min f(x), x G R n 

subject to 

Ax = b, Ae R m ' n , m <n, rank(A) = m. 

Let x\ be a feasible initial point. If we consider an iteration of the form Xj+i = xi — aidi, 
then the sequence Xi consists of feasible points iff 

Adi = (22) 

The general solution of (22) can be written using the ABS formula (12) 

di = Hl +1 q (23) 

In (23) matrix H m+ i depends on parameters Hi, W and V and q <G R n can also be seen as 
a parameter. Hence the general iteration generating feasible points is 

x i+ i = Xi- aiH^ +1 q. (24) 

The search vector is a descent direction if d T \7f(x) = q T H m+ iV f (x) > 0. This condition can 
always be satisfied by a choice of q unless H m+ iS/ f(x) = 0, implying, from the structure of 
Null(H m+ i), that V/(x) = A A for some A, hence that x«+i is a KT point with A the vector 
of Lagrange multipliers. If Xj+i is not a KT point, then we can generate descent directions 
by taking 

q = QH m+ iVf(x) (25) 

with Q symmetric positive definite. We obtain therefore a large class of methods with four 
parameter matrices (Hi , W, V, Q) . 

Some well-known methods in the literature correspond to taking in (25) W = I and 
building the Abaffian as follows. 
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(a) The reduced gradient method of Wolfe. H m+1 is built via the implicit LU method. 

(b) The projection method of Rosen. H m+ \ is built via the Huang method. 

(c) The Goldfarb and Idnani method. H rn+ \ is built via a modification of Huang method 

where H\ is a symmetric positive definite approximation of the inverse Hessian of f(x). 

To deal with linear inequality constraints there are two approaches in literature. 

• The active set method. Here the set of equality constraints is augmented with some 
inequality constraints, whose selection varies in the course of the process till the final set 
of active constraints is determined. Adding or cancelling a single constraint corresponds 
to a rank-one correction to the matrix defining the active set. The corresponding change 
in the Abaffian can be performed in order two operations, see Zhang (1995c). For an 
efficient ABS reformulation of the well-known active set method of Goldfarb ande Idnani 
see Xia, Liu and Zhang (1992) and Xia, Zhang and Liu (1995). 

• The standard form approach. Here one uses slack variables to put the problem in the 
following equivalent standard form 

min fix), 

with the constraints 

Ax = b, x>0. 

If x\ satisfies the above constraints, then a sequence of feasible points is generated by the 
following iteration 

x i+ i = Xi - ai(3iH m+1 Vf(x) (26) 

where a« may be chosen by a line search along vector H m+ iV f(x), while j3{ > is chosen to 
avoid violation of the nonnegativity constraints. 

If f(x) is nonlinear, then H m+ \ can usually be determined once for all at the first iteration, 
since generally V/(x) changes from point to point, allowing the determination of a new search 
direction. However if f(x) = c T x is linear, in which case we obtain the LP problem, then to 
get a new search direction we have to change H m+ \. We already observed in section 7 that 
the simplex method in the ABS formulation is obtained by constructing the matrix H m+ \ 
via the implicit LX method and at each step modifying one of the unit vectors used to build 
the Abaffian. One can show, see Xia (1995b), that the Karmarkar method (to be found in 
a previous paper by Evtushenko (1974)), corresponds to Abaffians built via a variation of 
the Huang algorithm, where the initial matrix is H\ = Diag(xi) and is changed at every 
iteration (whether the update of the Abaffian can be performed in order two operations is 
still an open question). Xia (1995b) discusses several questions related to the class of interior 
feasible point methods obtained via such a generalization of Huang algorithm. 
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9. Other results and conclusion 



Here we recall very briefly a few other results of ABS methods applied to optimization. 

• The computation of inertia of KT matrices is important in several algorithms recently 
proposed for quadratic programming. It can be done efficiently using a variation of the 
Huang algorithm, see Zhang (1999). 

• The LC problem can be reformulated with advantage using the implicit LX algorithm, 
see Feng et al. (1997). 

• ABS methods provide a new computationally interesting approach to concave program- 
ming, see Xia and Zhang (1999). 

• Explicit solution of some classes of integer linear inequalities and LP problems can be 
provided via the ABS algorithm for Diophantine equations, see Esmaeili and Mahdavi- 
Amiri (2000). 

Future work on ABS methods will not only deal with theoretical problems where ABS meth- 
ods can provide a fresh approach (e.g. the eigenvalue problem or even combinatorial prob- 
lems), but will be concentrated on the development of software, in the framework of the 
ABSPACK project. 
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RESULTS ON DETERMINED LINEAR SYSTEMS 
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RESULTS ON OVERDETERMINED SYSTEMS 
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RESULTS ON UNDERDETERMINED LINEAR SYSTEMS 
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RESULTS ON KT SYSTEMS 
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APPENDIX 2 



Vespucci et al. (1992) have given a 205 entries bibliography of ABS methods up to the end of 
1991, Nicolai and Spedicato (1997) have updated that bibliography to just over 300 entries. 
The present bibliography, tentatively complete up to the end of 1999, contains 349 entries 
(not counting earlier versions of some papers). 
The following abbreviations are used: 



JOTA 

OMS 

ABSC1 



ABSC2 



ABSDAL 
QDMSIA 



Journal of Optimization Theory and Applications 
Optimization Methods and Software 

Proceedings of the First International Conference on ABS 

Algorithms, Luoyang, September 1991, 

University of Bergamo (Spedicato E. editor), 1992 

Proceedings of the Second International 

Conference on ABS Algorithms, 

Beijing, June 1995, Report DMSIA 19/97, 

University of Bergamo 

(Spedicato E., Li Z. and Luksan L. editors) 

Collection of reports on the ABS algorithms, 

Dalian University of Technology, 1991 

Quaderni del Dipartimento di Matematica, Statistica, 

Informatica ed Applicazioni, University of Bergamo 



Abaffy J. (1979), A linearis egyenletrendszerek altalanos megoldasanak egy 
modszerosztalya, Alkalmazott Matematikai Lapok 5, 223-240. 

Abaffy J. (1984), A generalization of the ABS class, Proceedings of the Second Confer- 
ence on Automata Languages and Mathematical Systems, 7-11, University of Economics, 
Budapest. 

Abaffy J. (1986), Some special cases of the ABS class for band matrices, Proceedings of 
the Fourth Conference on Automata Languages and Mathematical Systems, 9-15, University 
of Economics, Budapest. 

Abaffy J. (1987a), Error analysis and stability for some cases in the ABSg class, Report 
190, Numerical Optimisation Centre, Hatfield Polytechnic. 

Abaffy J. (1987b), The description of a package for sparse linear systems using six algo- 
rithms in the ABSg class, Report 192, Numerical Optimisation Centre, Hatfield Polytechnic. 

Abaffy J. (1987c), Preliminary test results with some algorithms of the ABSg Class, 
Report 193, Numerical Optimisation Centre, Hatfield Polytechnic. 

Abaffy J. (1987d), Optimal methods in Broyden's stability sense, Report 194, Numerical 
Optimisation Centre, Hatfield Polytechnic. 
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Abaffy J. (1988a), A superlinear convergency theorem in the ABSg class for nonlinear 
algebraic equations, JOTA 59, 39-43. 

Abaffy J. (1988b), Equivalence of a generalization of Slobodans algorithm with a subclass 
of the generalized ABS algorithm for linear systems, QDMSIA 1/88. 

Abaffy J. (1988c), Derivation of Lanczos type methods in the generalized ABS algorithm 
for linear systems, QDMSIA 2/88. 

Abaffy J. (1988d), Formulation of the conjugate gradient type algorithm of Fridman for 
indefinite symmetric systems in the subclass of the generalized symmetric ABS algorithms, 
QDMSIA 3/88. 

Abaffy J. (1988e), Reorthogonalized methods in the ABSg class, QDMSIA 17/88. 

Abaffy J. (1991), ABS algorithms for sparse linear systems, in Computer Algorithms for 
Solving Linear Algebraic Equations: the State of the Art, NATO ASI Series, Vol. F77, 
111-132, Springer- Verlag (Spedicato E. editor). 

Abaffy J. (1995), A new method for computing eigenvalues, ABSC2. Abaffy J. and 
Bertocchi M. (1990), Vectorized ABS conjugate direction algorithms on the IBM 3090 VF, 
Pure Mathematics and Applications, Ser.B, 1, 189-202. 

Abaffy J. and Bertocchi M. (1991), A vectorized FORTRAN code of the implicit conjugate 
direction algorithm of the ABS class, Manual 1/57, IAC, Rome. 

Abaffy J., Bertocchi M. and Torriero A. (1991a), Characterization ofM-matrices via ABS 
methods and its application to input-output analysis, QDMSIA 11/91. 

Abaffy J., Bertocchi M. and Torriero A. (1991b), Una caratterizzazione delle M-matrici 
attraverso la classe ABS: una applicazione al modello input-output, Atti del XV Convegno 
dell'Associazione per la Matematica Applicata alle Scienze Economiche e Sociali, Grado, 
September 1991. 

Abaffy J., Bertocchi M. and Torriero A. (1992a), Criteria for transforming a Z-matrix 
into an M-matrix, OMS 1, 183-196. Also as QDMSIA 12/91. 

Abaffy J., Bertocchi M. and Torriero A. (1992b), Analysis of a subclass of nonsingular 
M-matrices and its economic applications, QDMSIA 22/92. 

Abaffy J., Bertocchi M. and Torriero A. (1992c), Un algoritmo per la determinazione 
dell' autovalore di Frobenius di matrici non negative, Atti XVI Convegno AMASES, Treviso, 
763-767. 

Abaffy J., Bertocchi M. and Torriero A. (1992d), Perturbations of M-matrices via ABS 
methods and their applications to input-output analysis, QDMSIA 23/92. Abaffy J., Bertoc- 
chi M. and Torriero A. (1992e), New algorithms for computing the Frobenius eigenvalue, 
QDMSIA 25/92. 

Abaffy J., Bertocchi M. and Torriero A. (1994), Sul calcolo delV autovalore di Frobenius 
di matrici non negative, Scritti in onore di Giovanni Melzi, Vita e Pensiero, Milano, 9-26 
(Manara C.F. and Faliva M. editors). 

Abaffy J., Broyden C.G. and Spedicato E. (1982), A class of direct methods for linear 
systems IF. generalization, non-singular representation and other tales, Report SOFMAT 
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21/82, IAC, Rome. 

Abaffy J., Broyden C.G. and Spedicato E. (1983), Numerical performance of the pseu- 
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Abaffy J., Galantai A. and Spedicato E. (1989), Forward error analysis in ABS methods, 
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translation with update appendix published in 1996 by Mir, Moscow). 

Abaffy J. and Spedicato E. (1989b), On the use of the ABS algorithm for some linear 
programming problems, Preprint, University of Bergamo. 

Abaffy J. and Spedicato E. (1990), A class of scaled direct methods for linear systems, 
Annals of the Institute of Statistical Mathematics, 42, 187-201. 

Adib M, Mahdavi Amiri N. and Spedicato E. (1999), Broyden's method as and ABS 
algorithm, Publ. Univ. Miskolc, Series D, 40, 3-13. 

Belyankov A.Y., (1990), Exhaustion methods in numerical linear algebra, Proceedings on 
International Conference on Optimization Software, Birkhauser 

Bertocchi M. (1989), Numerical experiments with ABS algorithms for linear systems on 
a parallel machine, JOTA 60, 375-392. 

Bertocchi M. (1991), A generalization of the orthogonally and optimally stable scaled 
ABS algorithms for linear systems, Pure Mathematics and Applications, Serie B, 2, 251-266. 

Bertocchi M. and Spedicato E. (1988), Vectorizing the modified Huang algorithm of the 
ABS class on the IBM 3090 VF, QDMSIA 23/88. Also in Proceedings of the International 
Meeting on Parallel Computing, Verona, Italy, 1988. 

Bertocchi M. and Spedicato E. (1989a), Performance of the implicit Gauss- Cholesky 
algorithm of the ABS class on the IBM 3090 VF, Proceedings of the 10-th Symposium on 
Algorithms, Strbske Pleso, 1989, 30-34. 

Bertocchi M. and Spedicato E. (1989b), Computational performance on the IBM 3090 
VF of the modified Huang and the implicit Gauss-Cholesky algorithms versus the Gaussian 
solver in the ESSL library on ill-conditioned problems, Proceedings of the 10-th Symposium 
on Algorithms, Strbske Pleso, April 17-21, 1989, 22-29. Also as QDMSIA 14/88. 

Bertocchi M. and Spedicato E. (1989c), Vectorizing the implicit Gauss-Cholesky algorithm 
of the ABS class on the IBM 3090 VF, Report 1/37, IAC, Rome. Also as QDMSIA 22/88. 

Bertocchi M. and Spedicato E. (1990a), Block ABS algorithms for dense linear systems 
in a vector processor environment, Proceedings of the Conference on Supercomputing Tools 
for Science and Engineering, Pisa, December 1989, Franco Angeli, Milano, 39-46 (Laforenza 
D. and Perego R. editors). Also as QDMSIA 2/90. 

Bertocchi M. and Spedicato E. (1990b), A vectorized FORTRAN code of the implicit 
Gauss-Cholesky algorithm of the ABS class, Manual 1/34, IAC, Rome. 

Bertocchi M., Brandolini L. and Spedicato E. (1990a), A vectorized FORTRAN code of 
the modified Huang algorithm of the ABS class, Manual 1/35, IAC, Rome. 

Bertocchi M., Brandolini L. and Spedicato E. (1990b), Vectorizing the modified Huang 
algorithm of the ABS class on the IBM 3090VF, Report 1/36, IAC, Rome. 

Bertocchi M., Spedicato E. and Vespucci M.T. (1991), A vectorized FORTRAN code of 
the implicit QR algorithm of the ABS class, Manual 1/56, IAC, Rome. 



22 



Bertocchi M., Spedicato E. and Vespucci M.T. (1992), Vectorizing the implicit QR al- 
gorithm of the ABS class on the IBM 3090 VF, Proceedings of the Meeting on Parallel 
Computing: Problems, Methods and Applications, Elsevier, 99-108 (Messina P. and Murli 
A. editors). 

Bodon E. (1987), Algoritmi ABS per minimi quadrati ed equazioni non lineari, Ph.D. 
Dissertation, University of Genova. 

Bodon E. (1989a), A code for linear full rank least squares problems based upon an 
implicit QR ABS algorithm, QDMSIA 14/89. 

Bodon E. (1989b), Globally optimally conditioned update in the ABS class for linear 
systems, QDMSIA 17/89. 

Bodon E. (1989c), Biconjugate algorithms in the ABS class I: alternative formulation and 
theoretical properties, QDMSIA 3/89. 

Bodon E. (1990a), A code for solving determined or underdetermined full or deficient rank 
linear systems based upon the optimally conditioned ABS algorithm, Preprint, University of 
Bergamo. 

Bodon E. (1990b), Numerical behaviour of the implicit QR algorithm of the ABS class 
for nonlinear systems, Preprint, University of Bergamo. 

Bodon E. (1992a), A code of the Hestenes-Stiefel type ABS algorithm, QDMSIA 16/92. 

Bodon E. (1992b), Numerical experiments with Gauss-ABS algorithms on tridiagonal 
systems of linear equations, QDMSIA 31/92. 

Bodon E. (1993a), Numerical results on the ABS algorithms for linear systems of equa- 
tions, QDMSIA 9/93. 

Bodon E. (1993b), Numerical experiments with ABS algorithms on upper banded systems 
of linear equations, Report TR/PA93/14, CERFACS, Toulouse. Also as QDMSIA 17/92. 

Bodon E. (1993c), Numerical experiments with ABS algorithms on banded systems of 
linear equations, Report TR/PA/93/13, CERFACS, Toulouse. Also as QDMSIA 18/92. 

Bodon E. (1993d), Computational performance of the implicit QR algorithm for linear 
least squares, QDMSIA 17/93. 

Bodon E. (1993e), Two codes of block ABS algorithms for linear systems of equations, 
Preprint, University of Bergamo. 

Bodon E. (1995a), Some experiments with the ABS-GMRES method, Preprint, University 
of Bergamo. 

Bodon E. (1995b), Two codes of the ABS-GMRES method, Preprint, University of Berg- 
amo. 

Bodon E. and Spedicato E. (1990a), The scaled ABS algorithm and its application to 
linear least squares: some numerical results, QDMSIA 3/90. 

Bodon E. and Spedicato E. (1990b), Numerical evaluation of the implicit LU, LQ and 
QU algorithms in the ABS class, QDMSIA 28/90. 

Bodon E. and Spedicato E. (1991), Factorized ABS algorithms for linear systems: deriva- 
tion and numerical results, QDMSIA 14/91. 



23 



Bodon E. and Spedicato E. (1992), On some STOD-ABS algorithms for large linear 
systems, QDMSIA 15/92. 

Broyden C.G. (1985), On the numerical stability of Huang's and related methods, JOTA 
47, 401-412. 

Broyden C.G. (1991), On the numerical stability of Huang's update, Calcolo, 28 (1991), 
303-311. Also as QDMSIA 18/89. 

Ceribelli C. (1989), Implementazione dell' algoritmo di Huang per sistemi lineari a bande, 
Ms. Dissertation, University of Bergamo. 

Chen F. and Xia Z. (1999), Roundoff error analysis of the ABS algorithm, QDMSIA 10/99 

Chen Z. and Deng N. (1992), A general algorithm for underdetermined linear systems, 
ABSC1. 

Chen Z. and Deng N. (1993), An algorithm with n-step quadratic convergency for non- 
linear systems, Numerical Mathematics, a Journal of Chinese Universities 15, 32-39 

Chen Z. and Deng N. (1995), A class of algorithms for quadratic programming via the 
ABS approach, ABSC2. 

Chen Z. and Deng N. (1997), Some algorithms for the convex quadratic programming 
problem via the ABS approach, OMS 8, 157-170 

Chen Z., Deng N. and Spedicato E. (1999), Truncated difference ABS-type algorithm for 
nonlinear systems of equations, OMS 10, 1-12 

Chen Z. and Spedicato E. (1993), Computing the null space of a matrix via the ABS 
approach, Preprint, University of Bergamo. 

Chen Z. and Xue Y. (1992), Efficient ABS-type algorithms for solving linear systems of 
equations, Preprint, Polytechnical University, Beijing. 

Deng N. (1998). A new nonlinear ABS method with convergence and efficiency analysis, 
QDMSIA 3/98. 

Deng N. and Chen Z. (1990), Truncated nonlinear ABS algorithm and its convergence 
property, Computing 45, 169-173. 

Deng N. and Chen Z. (1991), A unihed convergence theorem for nonlinear ABS-type 
algorithm, Preprint, University of Bergamo. 

Deng N. and Chen Z. (1998), A new nonlinear ABS-type algorithm and its efficiency 
analysis, OMS 10, 71-86 

Deng N. and Spedicato E. (1988), Optimal conditioning parameter selection in the ABS 
class through a rank two update formulation, QDMSIA 18/88. 

Deng N. and Xiao Y. (1989), The minimum norm correction class for linear system of 
equations, Preprint, Department of Applied Mathematics, Beijing Polytechnic University. 

Deng N. and Wang L. (1993), A modified ABS-type algorithm depending on nonlinearity, 
Journal of Beijing Agricultural Engineering University 13, 26-32. 

Deng N. and Zhu M. (1992a), Further results on the local convergence of the nonlinear 
ABS algorithm, OMS 1, 211-222. Also as QDMSIA 15/91. 



24 



Deng N. and Zhu M. (1992b), The local convergence of the nonlinear Voevodin-ABS 
algorithm, OMS 1, 223-231. Also as QDMSIA 7/91. 

Deng N., Li Z. and Spedicato E. (1996), On sparse Quasi-Newton quasi-diagonally dom- 
inant updates, QDMSIA 1/96. 

Deng N., Spedicato E. and Vespucci M.T. (1991), Experiments with the ABS implicit 
Gauss- Cholesky algorithm on nested dissection matrices, Report 1/69, IAC, Rome. 

Deng N., Spedicato E. and Zhu M. (1994), A local convergence theorem for nonlinear ABS 
algorithms, Computational and Applied Mathematics 13, 49-59. Also as QDMSIA 9/91. 

Deng N., Xiao Y. and Zhou F. (1992), Variationally derived ABS methods, OMS 1, 197- 
210. 

Dixon L.C.W. (1996), On the numerical stability of the implicit LU method and related 
methods, QDMSIA 18/97. 

Egervary E. (1953), On a property of the projector matrices and its application to canon- 
ical representation of matrix functions, Acta Scientiarum Mathematicarum 15, 1-6, 1953. 

Egervary E. (1957), Uber eine Methode zur numerische Losung der Poissonischen Differ- 
entialgleichung fur beliebige Gebiete, Acta Mathematica Academiae Scientiae Hungaricae 9, 
341-361. 

Egervary E. (1960), On rank-diminishing operations and their application to the solution 
of linear equations, ZAMP 9, 377-386. 

Esmaeili H. (1998), Application of the ABS algorithm to Diophantine equations and other 
tales, QDMSIA 4/98. 

Esmaeili H., N. Amiri-Mahdavi and Spedicato E. (1999), Solution of Diophantine systems 
via the ABS methods, QDMSIA 29/99 

Eremin Y. (1991), A block generalization of the ABS algorithm, Preprint, Institute of 
Numerical Analysis, Russian Academy of Science, Moscow. 

Faridpour H. (1994), A class of scaled direct methods for solving a system of linear 
equations, Ms. Dissertation, Tabriz University. 

Fegyverneki S. (1995), Some new estimators for robust nonlinear regression, Publications 
of University of Miskolc, Series D, 36, 23-32. 

Feng E. (1994), Implementing and testing an ABS algorithm for linearly constrained 
optimization, QDSMIA 17/94. 

Feng E. (1995), The FORTRAN code LCABSM of an ABS algorithm for linearly con- 
strained optimization, QDSMIA 1/95. 

Feng E. and Wang X. (1991), Application of the ABS algorithm to solve the mathematical 
simulation of geological history, ABSDAL. 

Feng E., Shi G. and Zhang L. (1994), The ABS formulation of weighted generalized 
inverse, Northeast Operations Research, 8, 39-41. 

Feng E., Wang X. and Wang X.L. (1995), The application of ABS algorithm in the dual 
simplex method, ABSC2. 



25 



Feng E., Wang X. and Wang X.L. (1997), On the application of the ABS algorithm to 
linear programming and linear complementarity, OMS 8, 133-142 

Fodor S. (1999a), A class of ABS methods for Diophantine systems of equations, part I, 
preprint, University of Economics, Budapest 

Fodor S. (1999b), A class of ABS methods for Diophantine systems of equations, part II, 
preprint, University of Economics, Budapest 

Fragnelli V. and Resta G. (1985), An ABS method for solving suitably structured linear 
systems is supported by parallel architectures, Proceedings of the First European Workshop 
on Parallel Processing Techniques for Simulation, Manchester, Plenum Press, 234-247. 

Galantai A. (1989), A new convergence theorem for nonlinear ABS methods, QDMSIA 
12/89. 

Galantai A. (1990a), Block ABS methods for nonlinear systems of algebraic equations, 
in Computational Solution of Nonlinear Systems of Equations, SIAM Lectures in Applied 
Mathematics 26, 181-187 (Allgower E.L. and Georg K. editors). 

Galantai A. (1990b), A global convergence theorem for nonlinear ABS methods, QDMSIA 
9/90. 

Galantai A. (1991a), Convergence theorems for the nonlinear ABS methods, Colloquia 
Mathematica Societatis Janos Bolyai 59, Numerical Methods, Miskolc 1990, North-Holland, 
65-79 (Stoyan G. editor). 

Galantai A. (1991b), Analysis of error propagation in the ABS class, Annals of the Insti- 
tute of Statistical Mathematics 43, 597-603. Also as QDMSIA 1//87. 

Galantai A. (1993a), Generalized implicit LU algorithms in the class of ABS methods for 
linear and nonlinear systems of algebraic equations, QDMSIA 5/93. 

Galantai A. (1993b), ABS methods on large nonlinear systems with banded Jacobian 
structure, QDMSIA 14/93. 

Galantai A. (1993c), Testing of implicit LU ABS methods on large nonlinear systems with 
banded Jacobian, QDMSIA 19/93. 

Galantai A. (1994), ABS methods for the solution of linear and nonlinear systems of 
equations, Habilitation Dissertation, Technical University of Budapest. 

Galantai A. (1995a), Hybrid Quasi-Newton methods for nonlinear equations, ABSC2. 

Galantai A. (1995b), The global convergence of the ABS methods for a class of nonlinear 
problems, OMS 4, 283-295. Also as QDMSIA 9/92. 

Galantai A. (1995c), A fast method for solving systems of nonlinear algebraic equations, 
Publications of the University of Miskolc, Series C, 45, 85-92. 

Galantai A. (1997a), Hybrid Quasi-Newton ABS methods for nonlinear equations, QDM- 
SIA 11/97. 

Galantai A. (1997b), The local convergence of Quasi-Newton ABS methods, Publications 
of the University of Miskolc, Series D, 37, 31-38 

Galantai A. (1999a), Perturbation theory for full rank factorizations, QDMSIA 40/99 



26 



Galantai A. (1999b), Parallel ABS projection methods for linear and nonlinear systems 
with block arrohead structure, Computers and Mathematics with Applications 38, 11-17 

Galantai A. and Jeney A. (1993), Quasi-Newton ABS methods, Proceedings of MICROCAD- 
93 Conference, Miskolc, 63-68. 

Galantai A. and Jeney A. (1996), Quasi-Newton ABS methods for solving nonlinear al- 
gebraic systems of equations, JOTA 89, 561-573. 

Galantai A., Jeney A. and Spedicato E. (1993a), Testing of ABS-Huang methods on 
moderately large nonlinear systems, QDMSIA 6/93 

Galantai A., Jeney A. and Spedicato E. (1993b), A FORTRAN code of a Quasi-Newton 
ABS algorithm, QDMSIA 7/93. 

Ge R. (1997), A family of Broyden-ABS type methods for solving nonlinear systems, 
QDMSIA 1/97. 

Ge R. and Xia Z. (1996), An ABS algorithm for solving nonsingular systems, QDMSIA 
36/96. 

Gecan A. and Snider A.D. (1997), Rationale for the Camerini-Fratta-Maffioli modified 
subgradient deflection and its ABS formulation, ABSC2, 57-66 

Gredzhuk V.I. and Petrina D.Y. (1998), The ABS method and the method of averaging 
functional corrections, QDMSIA 8/98 

Giudici S. (1993), Applicazione del metodo ABS per la soluzione dei problem! di minimi 
quadrati nonlineari, Ms. Dissertation, University of Bergamo. 

Gu S. (1989), The simplified ABS algorithms, QDMSIA 5/89. 

Huang H.Y. (1975), A direct method for the general solution of a system of linear equa- 
tions, JOTA 16, 429-445. 

Huang Z. (1989), On the ABS methods for nonlinear problems, Ph.D. Dissertation, Fudan 
University, Shanghai. 

Huang Z. (1991a), Multi-step nonlinear ABS methods and their efficiency analysis, Com- 
puting 46, 143-153. 

Huang Z. (1991b), Modified ABS methods for nonlinear systems without evaluating Ja- 
cobian matrices, Numerical Mathematics, a Journal of Chinese Universities 13, 60-71. 

Huang Z. (1992a), Convergence analysis of the nonlinear block scaled ABS methods, 
JOTA 75, 331-344. Also as QDMSIA 14/90. 

Huang Z. (1992b), A family of discrete ABS type algorithms for solving systems of non- 
linear equations, Numerical Mathematics, a Journal of Chinese Universities 14, 130-143. 

Huang Z. (1992c), ABS methods for solving nonlinear least squares problems, Communi- 
cations on Applied Mathematics and Computation 6, 75-85. 

Huang Z. (1993a), A class of ABS methods for nonlinear least squares problems, Ricerca 
Operativa 26, 67-78. Also as QDMSIA 3/92. 

Huang Z. (1993b), On the convergence analysis of the nonlinear ABS methods, Chinese 
Annals of Mathematics 14 B, 213-224. 



27 



Huang Z. (1993c), Row update ABS methods for solving sparse nonlinear systems of 
equations, OMS 2, 297-309. 

Huang Z. (1993d), Restart row update ABS methods for solving systems of nonlinear 
equations, Computing 50, 229-239. Also as QDMSIA 13/92. 

Huang Z. (1994a), Row update ABS methods for systems of nonlinear equations, OMS 
3, 41-60. Also as QDMSIA 14/92. 

Huang Z. (1994b), Generalization of row update ABS methods for systems of nonlinear 
equations, QDMSIA 14/94. 

Huang Z. (1994c), A new method for solving nonlinear underdetermined systems, Com- 
putational and Applied Mathematics 1, 33-48. Also as QDMSIA 5/92. 

Huang Z. and Spedicato E. (1994d), Numerical experiments with Newton's, Brent's and 
Huang type methods for nonlinear systems, QDMSIA 1/94. 

Huang Z. and Spedicato E. (1994e), Numerical testing of Quasi-Newton and some related 
methods for nonlinear systems, QDMSIA 4/94. 

Huang Z. and Spedicato E. (1994f), Further numerical testing of Newton-type methods 
for nonlinear systems, QDMSIA 9/94. 

Huang Z. and Spedicato E. (1996), ABS formulation of the GMRES method, QDMSIA 
14/96. 

Jeney A. (1991), Numerical studies of ABS methods for solving nonlinear systems of 
equations, Alkalmazott Matematikai Lapok, 15, 353-364 (in Hungarian). 

Jeney A. (1991), Discretization of a subclass of ABS methods for solving nonlinear systems 
of equations, Alkalmazott Matematikai Lapok, 15, 365-379 (in Hungarian). Also as QDMSIA 
10/92 (in English). 

Jeney A. (1992a), The performance of ABS-methods for solving nonlinear systems of 
algebraic equations, QDMSIA 11/92. 

Jeney A. (1992b), Discrete ABS methods for nonlinear systems of equations, Proceedings 
of MICROCAD-92, Miskolc, 107-110 (in Hungarian). 

Jeney A. (1995), The convergence of discrete ABS methods and numerical investigations, 
Proceedings of MICROCAD-95, Miskolc, Section K, 22-25. 

Jeney A. (1996a), Calculation of deformations of fibre-reinforced hoses with system of 
nonlinear algebraic equations, Proceedings of MICROCAD-96, Section L, Miskolc, 25-30. 

Jeney A. (1996c), Numerical comparison of multistep Quasi-Newton ABS methods with 
the multistep Broyden method for solving systems of nonlinear equations, Publications of the 
University of Miskolc, Series D, 36, 47-54 

Jeney A. (1997a), The local convergence of Quasi-Newton ABS methods , Publications 
of the University of Miskolc, Series D, 37, 31-38. 

Jeney A. (1997b), Multistep discrete ABS methods for solving systems of nonlinear equa- 
tions, Publications of the University of Miskolc, Series D, 37, 39-46. 

Lai Y., Gao Z. and He G. (1995), An ABS-generalized projection algorithm for linearly 
constrained optimization problems, ABSC2. 



28 



Li C. (1997), On the preconditioning properties of an ABS-CG type algorithm for solving 
large linear systems, International Conference on Computational Mathematics, Bangkock 

Li C. (1999), si On a new algorithm of ABS-CG type for solving nonlinear systems, 
Numerical Methods and Computer Applications 20 

Li S. (1997), The ABS methods for rank one perturbed linear system of equation, ABSC2, 
45-52 

Li S. (1995), On a variational characterization of the ABS algorithm, preprint. 
Li S. (1995), On the use of an ABS method in Karmarkar algorithm, preprint. 
Li S. and Zhang L. (1991), Some methods for base matrix Null(A) and numerical exper- 
iments, ABSDAL. 

Li Z.F. (1996), Restarted and truncated versions of ABS methods for large linear systems: 
a basic framework, QDMSIA 8/96. 

Li Z.F. (1996), Modifying the Huang algorithm for parallel computation of solutions to 
large and unstructured linear systems, QDMSIA 25/96. 

Li Z.F. (1996), Truncated block ABS methods for large linear systems, QDMSIA 26/96. 

Li Z.F. (1996), Applications of ABS methods to infinite linear systems, QDMSIA 27/96. 

Li Z.F. (1996), On the convergence of truncated ABS algorithms for large linear systems, 
QDMSIA 29/36. 

Li Z.F. (1996), Towards ABS methods for ill-posed equality constrained linear least 
squares problems, QDMSIA 31/96. 

Li Z.F., Deng N. and Wang L. (1996) ABS-type methods with ordering strategies for 
solving nonlinear equations, QDMSIA 12/96. 

Longhino M. (1991), Un algoritmo ABS di tipo gradient! coniugati per minimi quadrati 
lineari, Ms. Dissertation, University of Milan. 

Longhino M., Spedicato E. and Vespucci M.T. (1992), On the minimum residual conjugate 
gradient type ABS algorithm for linear least squares, QDMSIA 29/92. 

Luksan L. and Spedicato E. (1998) Variable metric methods for unconstrained optimiza- 
tion, QDMSIA 7/98 

Luo Y., Feng Z. and Liu J. (1997), On the conditions of the semilocal convergence for 
nonlinear ABS algorithms, ABSC2, 27-34 

Meng F., Wang X. and He Z. (1997), A theorem on W m being K, L admissible, ABSC2, 
53-56 

Mirnia K. (1996), Iterative rehnement in ABS methods, QDMSIA 32/96. 

Nicolai S. (1995), Metodi numerici per la soluzione delle equazioni KT, Ms. Dissertation, 
University of Bergamo. 

Nicolai S. and Spedicato E. (1997), A bibliography of the ABS methods, OMS 8, 171-183 

Nicolai S. and Spedicato E. (1996), A bibliography on the ABS methods, QDMSIA 3/96. 

Oprandi M. (1987), Metodi ABS per sistemi algebrici lineari, Ms. Dissertation, University 
of Bergamo. 



29 



PhuaK.H. (1986), Solving sparse linear systems by ABS methods, Report 4/86, Computer 
Science Department, National University, Singapore. 

Phua K.H. (1988), Solving sparse linear systems by an ABS method that corresponds to 
LU decomposition, BIT 28, 709-718. 

Ren Z. (1991), Solving optimization problems by combining the ABS methods and mono- 
tonic analysis, ABSDAL. 

Ronto N.I., Tuzson A. (1994) Construction of periodic solutions of differential equations 
with impulse effect, Publ. Math. Debrecen, 44, 335-357. 

Shahmorad S. (1994), Solution of linear least squares via the ABS algorithm, Ms. Disser- 
tation, Tabriz University. 

Shan R. and Liu W. (1997), A hnite terminating conjugate direction method for linearly 
constrained optimization, ABSC2, 35-44 

Shi G. (1992), An ABS algorithm for generating nonnegative solutions to linear systems, 
ABSC1, 54-57. Also in ABSDAL. 

Shi G. and Liu X. (1992), An ABS algorithm for the concave-minimum problem over a 
convex polyhedroid, ABSC1, 58-64. Also in ABSDAL. 

Spedicato E. (1983), On the symmetric algorithm in the ABS class of direct methods for 
linear systems, QDMSIA 13/83. 

Spedicato E. (1985), On the solution of linear least squares through the ABS class for 
linear systems, Atti Giornate AIRO, September 1985, 89-98. 

Spedicato E. (1986), Proprieta di convergenza degli algoritmi ABS per sistemi algebrici 
nonlineari, Ricerca Operativa e Informatica, Milano, Franco Angeli, 57-65 (Bielli M. editor). 

Spedicato E. (1987a), A subclass of variationally derived algorithms in the ABS class for 
linear systems, Report 202, Mathematische Institute, University of Wiirzburg. 

Spedicato E. (1987b), Optimal conditioning parameter selection in the ABS class for 
linear systems, Report 203, Mathematische Institute, University of Wiirzburg. 

Spedicato E. (1987c), ABS algorithms for sequential and parallel solution of linear and 
nonlinear algebraic equations, QDMSIA 12/87. 

Spedicato E. (1987d), Variationally derived algorithms in the ABS class for linear systems, 
Calcolo, 24, 241-246. 

Spedicato E. (1988), Optimal conditioning parameter selection in the ABS class for linear 
algebraic systems, Matematica Aplicada e Computacional, 7, 187-199. 

Spedicato E. (1989), A class of algorithms for sequential and parallel solution of algebraic 
linear and nonlinear systems, Advances in Chinese Mathematics 18, 55-61 (in Chinese). Also 
in Proceedings of the 7th Mathematical Programming Symposium, Nagoya, 1986 (in English). 

Spedicato E. (1991a), ABS algorithms for general linear systems, in Computer Algorithms 
for Solving Linear Algebraic Equations: the State of the Art, NATO ASI Series, Vol. F77, 
93-110, Springer- Verlag (Spedicato E. editor). 

Spedicato E. (1991c), ABS algorithms for general linear systems, QDMSIA 8/91. 



30 



Spedicato E. (1992a), A class of sparse symmetric quasi-Newton updates, Ricerca Oper- 
ativa 22, 63-70, 1992. Also as QDMSIA 1/92. 

Spedicato E. (1992b), Special Issue: The First International Conference on ABS Methods, 
OMS 1, 169-281 (editor). 

Spedicato E. (1992c), Proceedings of the First International Conference on ABS algo- 
rithms, University of Bergamo (editor). 

Spedicato E. (1992d), A review of ABS algorithms for general linear systems, OMS 1, 
171-182. 

Spedicato E. (1993), Ten years of ABS methods: a review of theoretical results and 
computational achievements, Surveys on Mathematics for Industry, 3, 217-232. 

Spedicato E. (1995a), ABS algorithms for linear systems and linear least squares: the- 
oretical results and computational performance, Scientia Iranica 1, 289-303, 1995. Also as 
QDMSIA 2/94. 

Spedicato E. (1995b), ABS algorithms from Luoyang to Beijing, QDMSIA 12/95. 

Spedicato E. (1996), New ABS linear solvers and applications to optimization, PhD Dis- 
sertation, Dalian University of Technology. 

Spedicato. (1997), ABS algorithms from Luoyang to Beijing, OMS 8, 87-97 

Spedicato E. (1998), De aliqua method! Huangianae proprietate, Vox Latina 347-351, also 
as QDMSIA 4/91. 

Spedicato E. (1998), ABS algorithms for linear equations and linear least squares: a 
review, QDMSIA 15/98. 

Spedicato E. (1999), ABS algorithms for linear equations and optimization: a review of 
main results, Proceedings ISC99, Beirut, March 1999 (I. Moghrabi and S. Kabbani editors), 
LAU University, 11-21 

Spedicato E. and Abaffy J. (1982), A class of direct methods for linear systems I: basic 
properties, Report 82/4, IAMI, Milan. 

Spedicato E. and Abaffy J. (1987) On the use of the ABS algorithms for some linear 
programming problems, Preprint, University of Bergamo. 

Spedicato E. and Bodon E. (1987a), Numerical experiments in the ABS class for nonlinear 
systems of algebraic equations, Report SOFMAT 1/87, IAC, Rome. 

Spedicato E. and Bodon E. (1987b), Application of the ABS algorithm to linear least 
squares, QDMSIA 13/87. 

Spedicato E. and Bodon E. (1989a), Solving linear least squares by orthogonal factor- 
ization and pseudoinverse computation via the modified Huang algorithm in the ABS class, 
Computing, 42, 195-205. Also as QDMSIA 9/88. 

Spedicato E. and Bodon E. (1989b), Biconjugate algorithms in the ABS class II: numerical 
evaluation, QDMSIA 4/89. 

Spedicato E. and Bodon E. (1990), Numerical experiments of the implicit QR algorithm 
of the ABS class for nonlinear least squares, Preprint, University of Bergamo. 



31 



Spedicato E. and Bodon E. (1991a), Numerical performance of conjugate gradient type 
ABS methods on ill-conditioned problems, QDMSIA 13/91. 

Spedicato E. and Bodon E. (1991b), Numerical behaviour of the implicit QR algorithm 
in the ABS class for nonlinear systems, QDMSIA 16/91. 

Spedicato E. and Bodon E. (1992), Numerical behaviour of the implicit QR algorithm in 
the ABS class for linear least squares, Ricerca Operativa 22, 43-55. Also as QDMSIA 13/89. 

Spedicato E. and Bodon E. (1993), Solution of linear least squares via the ABS algorithm, 
Mathematical Programming 58, 111-136. Also as QDMSIA 11/89. 

Spedicato E., Chen Z. and Bodon E. (1996), ABS methods for KT equations, in Nonlinear 
Optimization and Applications, Plenum Press (Di Pillo G. and Giannessi F. editors). Also 
as QDMSIA 13/95. 

Spedicato E., Chen Z. and Deng N. (1993), The truncated difference ABS-type algorithm 
for nonlinear systems of equations, QDMSIA 21/93. 

Spedicato E., Chen Z. and Deng N. (1994), A class of difference ABS-type algorithms for 
a nonlinear system of equations, Numerical Linear Algebra with Applications 13, 313-329. 
Also as QDMSIA 6/91. 

Spedicato E., Chen Z. and Deng N. (1996), A new nonlinear ABS-type algorithm and its 
convergence property, QDMSIA 33/96. 

Spedicato E., Gao Z. and Yu H. (1990), The scaled ABS versions of some conjugate 
gradient methods for nonsymmetric linear systems, QDMSIA 19/90. 

Spedicato E. and Huang Z. (1992), An orthogonally scaled ABS method for nonlinear 
least squares problems, OMS 1, 233-242. Also as QDMSIA 4/92. 

Spedicato E. and Huang Z. (1995), Optimally stable ABS methods for nonlinear under- 
determined systems, OMS 5, 17-26. Also as QDMSIA 20/92. 

Spedicato E. and Huang Z. (1996), Numerical experience with Newton-like methods for 
nonlinear algebraic equations, QDMSIA 96/15. 

Spedicato E. and Li Z. (1997), A class of sparse symmetric Quasi-Newton updates with 
least-change column-correction, QDMSIA 97/17. 

Spedicato E. and Li Z. (1998), Restarted and truncated versions of ABS methods for 
large linear systems: a basic framework, in Proceedings of the First Workshop on Large Scale 
Scientific Computations, Varna, June 7-11, 1997 

Spedicato E. and Li Z. (1998), Towards ABS methods for large-scale sparse linear systems, 
QDMSIA 98/1. 

Spedicato E., Li Z. and Luksan L. (1997), Proceedings of the second international con- 
ference on ABS algorithms. Beijing, June 1995, QDMSIA 97/19. 

Spedicato E. and Tuma M. (1993), Solving sparse unsymmetric linear systems by implicit 
Gauss algorithm: stability, QDMSIA 4/93. 

Spedicato E. and Vespucci M.T. (1987), Computational performance of the implicit Gram- 
Schmidt (Huang) algorithm for linear algebraic systems, QDMSIA 4/87. 



32 



Spedicato E. and Vespucci M.T. (1992), A bibliography on the ABS methods, QDMSIA 
19/92. 

Spedicato E. and Vespucci M.T. (1993), Variations on the Gram-Schmidt and the Huang 
algorithms for linear systems: a numerical study, Applications of Mathematics 2, 81-100. 
Also as QDMSIA 21/89. 

Spedicato E. and Xia Z. (1990a), Application of ABS algorithms to constrained optimiza- 
tion I: linear equality constraints, QDMSIA 7/90. 

Spedicato E. and Xia Z. (1990b), An approach for polyhedral approximations to subdif- 
ferentials of hnite convex functions, QDMSIA 17/90. 

Spedicato E. and Xia Z. (1992a), Finding general solutions of the Quasi-Newton equation 
via the ABS approach, OMS 1, 243-252. Also as QDMSIA 10/91. 

Spedicato E. and Xia Z. (1992b), A class of descent algorithms for nonlinear function 
minimization with linear equality constraints, OMS 1, 265-272. 

Spedicato E. and Xia Z. (1992c), On some ABS methods for the computation of the 
pseudoinverse, Ricerca Operativa 22, 35-41. Also as QDMSIA 24/90. 

Spedicato E. and Xia Z. (1994), ABS algorithms for nonlinear optimization, in Algorithms 
for Continuous Optimization: the State of the Art, Kluwer, 333-356 (Spedicato E. editor). 

Spedicato E. and Xia Z. (1998), ABS algorithms for optimization: a review, QDMSIA 
16/98. 

Spedicato E. and Xia Z. (1999), ABS formulation of Fletcher implicit LU method for 
factorizing a matrix and solving linear equations and LP problems, QDMSIA 1/99 

Spedicato E. and Yang Z. (1990), Optimally conditioned scaled ABS algorithms for linear 
systems, JOTA 67, 141-150. Also as QDMSIA 15/89 

Spedicato E. and Zhao J. (1992), On sparse symmetric and positive dehnite quasi-Newton 
updates, QDMSIA 2/92. 

Spedicato E. and Zhao J. (1993), Explicit general solution of the Quasi-Newton equation 
with sparsity and symmetry, OMS 2, 311-319. Also as QDMSIA 11/93. 

Spedicato E. and Zhu M. (1994), On the generalized implicit LU algorithm of the ABS 
class for linear systems, QDMSIA 3/94. 

Spedicato E. and Zhu M. (1996), Reformulation of the ABS algorithm via full rank Abaf- 
hans, QDMSIA 4/96. 

Spedicato E. and Zhu M. (1996), Efficient solution of linear systems after rank-one cor- 
rection via the ABS algorithm with full rank AbafRans, QDMSIA 5/96. 

Spedicato E. and Zhu M. (1999), A generalization of the implicit LU algorithm to an 
arbitrary initial matrix, Numerical Algorithms 20, 343-351 

Spedicato E., Esmaeili H. and Xia Z. (1999), A review of ABS algorithms for linear 
real and Diophantine equations and optimization, Proceedings XVI CEDYA and VI CMA 
Conferences, Las Palmas de Gran Canaria, September 1999 (Montenegro R., Montero G. and 
Winter G. editors), CEANI, 11-23 



33 



Spedicato E., Xia Z. and Gao Z. (1990), The basic block ABS method I: formulation and 
theoretical properties, QDMSIA 15/90. 

Spedicato E., Xia Z. and Zhang L. (1997), The implicit LX method of the ABS class, 
OMS 8, 99-110, also as QDMSIA 2/96. 

Spedicato E., Xia Z., Zhang L. and Mirnia K. (1996), ABS algorithms for linear equations 
and applications to optimization, QDMSIA 30/96. 

Spedicato E., Xia Z. and Zhang L. (1999), ABS algorithms for solving linearly constrained 
optimization problems via the active set strategy, QDMSIA 16/99 

Spedicato E., Zhang L. and Xia Z. (1999), Solving linear matrix equations via the ABS 
method and applications to Quasi-Newton equations, QDMSIA 12/99 

Tan Z. (1991), A new algorithm with finite convergence to solve a class of quadratic 
programming constructed by the ABS methods, ABSDAL. 

Tan Z. and Chen Y. (1991), An improvement on sub optimization methods for convex 
programming constraints, ABSDAL. 

Tan Z. and Gao Y. (1991), A quasi-gradient projection method for nonlinear constrained 
systems, ABSDAL. 

Tan Z. and Gao Y. (1992), An improved gradient projection method for nonlinear con- 
strained systems, ABSC1, 65-77. 

Tuma M. (1990), Utilization of matrix sparsity in ABS methods for solving systems of 
linear equations, QDMSIA 21/90. 

Tuma M. (1993), Solving sparse unsymmetric sets of linear equations based on implicit 
Gauss projections, Report 556, Institute of Computer Science, Academy of Sciences, Prague. 

Tuma M. (1996), On orthogonalization procedures in some iterative algorithms for systems 
of linear equations, QDMSIA 8/96. 

Vespucci M.T. (1990), On the implicit Gauss-Cholesky algorithm of the ABS class, Con- 
tributed Papers, NATO ASI on Computer Algorithms for Solving Linear Algebraic Equations: 
the State of the Art, II Ciocco, 1990, Cooperativa Studium Bergomense, 186-205 (Spedicato 
E. editor). 

Vespucci M.T. (1991a), The use of the ABS algorithms in truncated Newton methods for 
linear optimization, Colloquia Mathematica Societatis Janos Bolyai 59, Numerical Methods, 
81-91 (Stoyan G. editor). 

Vespucci M.T. (1991b), An ABS implicit Gauss-Cholesky code for ND matrices, Manual 
1/68, Progetto Finalizzato Sistemi Informatici e Calcolo Parallelo, IAC, Roma. 

Vespucci M.T. (1991c), Truncated Newton methods based on the ABS class, Ph.D. Dis- 
sertation, Hatfield Polytechnic. 

Vespucci M.T. (1992), The ABS preconditioned conjugate gradient algorithm, QDMSIA 
6/92. Also in ABSC1, 78-94. 

Vespucci M.T., Yang Z., Feng E., Yu H. and Wang X. (1992), A bibliography on the ABS 
methods, OMS 1, 273-281. 



34 



Vitali A. (1989), Algoritmi vettoriali e paralleli della classe ABS per la soluzione di sistemi 
lineari, Ms. Dissertation , University of Milan. 

Xia Z. (1990a), Application of ABS algorithms to constrained optimization II: linear 
inequality constraints, QDMSIA 29/90. 

Xia Z. (1990b), An algorithm for minimizing a class of quasidifferentiable functions, QDM- 
SIA 12/90. 

Xia Z. (1991a), Quadratic programming via the ABS algorithm I: equality constraints, 
QDMSIA 1/91. 

Xia Z. (1991b), The ABS class and quadratic programming II: general constraints, QDM- 
SIA 5/91. 

Xia Z. (1992a), Finding subgradients or descent directions of convex functions by external 
polyhedral approximation of subdifferentials, OMS 253-264. 

Xia Z. (1992b), Derivation of some linear transformations via the ABS algorithm, ABSC1, 
95-110. 

Xia Z. (1994), An efficient ABS algorithm for linearly constrained optimization, QDMSIA 
15/94. 

Xia Z. (1995a), ABS formulation and generalization of the reduced gradient method, 
QDMSIA 3/95. 

Xia Z. (1995b), ABS formulation and generalization of the interior point method, ABSC2. 

Xia Z. (1995c), ABS reformulation of some versions of the simplex method for linear 
programming, QDMSIA 10/95. 

Xia Z. (1998), Further ABS formulation of reduced gradient methods, QDMSIA 2/98. 

Xia Z. and Gao Z. (1990), The basic block ABS methods II: application to the Rutherford 
matrix equation, QDMSIA 16/90. 

Xia Z. and Liu Y. (1993), On ABS descent directions for constrained minimization with 
linear equality constraints, Journal of Dalian University of Technology 35, 255-263. 

Xia Z. and Zhang L. (1999), A modified outer approximation method for solving concave 
programming problems via the ABS algorithm, QDMSIA 21/99 

Xia Z., Liu Y. and Zhang L. (1992), Application of a representation of ABS updating 
matrices to linearly constrained optimization I, Northeast Operational Research 7, 1-9. 

Xia Z., Zhang L. and Liu Y. (1995) An ABS algorithm for minimizing a nonlinear function 
subject to linear inequalities, ABSC2. 

Xiu N. and Chen G. (1991), Optimality of the ABS matrices update formula, Preprint, 
Qufu Normal University. 

Xue Y. (1997), A practical method for solving linear and nonlinear least squares problems, 
ABSC2, 17-26 

Yang Z. (1989a), ABS algorithms for solving certain systems of indefinite equations, 
QDMSIA 6/89. 

Yang Z. (1989b), The properties of iterative point in ABS algorithms and the related 
generalized inverse matrix, QDMSIA 7/89. 



35 



Yang Z. (1989c), On the numerical stability of the Huang and the modified Huang algo- 
rithms and related topics, QDMSIA 8/89. 

Yang Z. (1989d), An ABS subclass with optimal conditioning parameter selection, QDM- 
SIA 10/89. 

Yang Z. and Zhu M. (1994), The practical ABS algorithm for large scale nested dissection 
linear system, QDMSIA 11/94. 

Wang X., Cai J. and Li X. (1995), On the application of the ABS algorithm in the simplex 
method, ABSC2. 

Wang X., Feng E. and Liu X. (1995), ABS reformulation of Lemke's method for the linear 
complementarity problem, ABSC2. 

Wang Z. (1997), ABS formulation of Bunch-Parlett decomposition, ABSC2, 9-16 

Zeng K. and Jiang G. (1991), An application of the ABS algorithms to solve nonlinear 
systems of equations, Preprint, South- Western Jiaotong University. 

Zeng L. (1989), An ABS algorithm for sparse linear systems, QDMSIA 9/89. 

Zeng L. (1991), The parallel scaled ABS class for sparse systems of linear equations, 
ABSDAL. 

Zeng L. (1997), A parallel ABS method for sparse nonlinear algebraic equations, ABSC2, 

1-8 

Zhang J. and Zhu M. (1995a), A symmetric ABS algorithm for some matrix matrix 
equation and its minimum properties, Report MA95-19, City University, Hong Kong. 

Zhang J. and Zhu M. (1995b), Solving a least change problem under the weak secant 
condition by an ABS approach, Report MA95-20, City University, Hong Kong 

Zhang J. and Zhu M. (1997), Least-change properties of ABS methods and their applica- 
tion in secant-type updates, OMS 8, 111-132 

Zhang L. (1992a), Application of the ABS algorithm to Quasi-Newton methods for linearly 
constrained optimization problems, ABSC1 119-130. Also in ABSDAL. 

Zhang L. (1992b), A method for finding a feasible point of inequalities, ABSC1 131-137. 
Also in ABSDAL. 

Zhang L. (1992c), A class of quasi-Newton formulae with ZZ T factorization, Proceedings 
of the First Youth Academic Meeting of Liaoning, 1992, 220-224, Northeast University Press. 

Zhang L. (1995a), An algorithm for the least Euclidean norm solution of a linear system 
of inequalities via the Huang ABS algorithm and the Goldfarb-Idnani strategy, QDMSIA 
2/95. 

Zhang L. (1995b), Relations between projection matrices and ABS updating matrices, 
QDMSIA 5/95. 

Zhang L. (1997), On the ABS algorithm with singular initial matrix and its application 
to linear programming, OMS 8 143-156 (also as QDMSIA 8/95). 

Zhang L. (1995d), Updating of AbafRan matrices under perturbation in W and A, QDM- 
SIA 15/95. 



36 



Zhang L. (1999), Computing inertias of KKT matrix and reduced Hessian via the ABS 
algorithm for applications to quadratic programming, QDMSIA 7/99 

Zhang L., Feng E. and Jiang Z. (1992), A class of projection methods for linearly con- 
strained optimization problems, ABSC1, 138-149. Also in ABSDAL. 

Zhang L., Feng E. and Xia Z. (1993), A theorem on the steepest descent direction for 
linearly constrained problem, Journal of Mathematical Research and Exposition 13, 358. 

Zhang L. and Xia Z. (1993), A modihed version of ABS algorithms, Journal of Mathe- 
matical Research and Exposition 13, 26. 

Zhang L. and Xia Z. (1994), Comparing ABS algorithm with Fletcher and GI methods, 
Northeast Operations Research, 8, 36-38. 

Zhang L. and Xia Z.H. (1995), Application of the implicit LX algorithm to the simplex 
method, QDMSIA 9/95. 

Zhang L., Xia Z. and Feng E. (1998), Introduction to ABS Methods in Optimization, 
Monograph, Dalian University of Technology Press 

Zhang L., Xia Z. and Wang X. (1994), A version of updating ABS matrices in optimiza- 
tion, Northeast Operations Research, 8, 51-54. 

Zhang J. and Zhu M. (1997), Least change properties of ABS methods and their applica- 
tion in secant-type updates, OMS 8, 111-132 

Zhao J. (1981), Huang's method for solution of consistent linear equations and its gener- 
alization, Journal of Numerical Mathematics of Chinese Universities, 3, 8-17. 

Zhao J. (1989), A class of direct methods for solving linear inequalities, Journal of Nu- 
merical Mathematics of Chinese Universities 3, 231-38. 

Zhao J. (1991), ABS algorithms for solving linear inequalities, QDMSIA 21/91. 

Zhao J. (1995), Anoteon the ABS algorithm for solving linear systems, Journal of Nanjing 
University 12, 32-36. 

Zheng H., Fei P. and Fan R. (1991), The parallel Huang algorithm of the ABS class on 
MIMD parallel system, Preprint, University of Wuhan. 

Zhu M. (1987), The implicit LL T algorithm for sparse nested dissection linear systems, 
Technical Report 196, Numerical Optimisation Centre, Hatfield Polytechnic. 

Zhu M. (1989a), The Q-quadratic convergence for a subclass of nonlinear ABS algorithms, 
QDMSIA 22/89. 

Zhu M. (1989b), A generalization of nonlinear ABS algorithms and its convergence prop- 
erties, Journal of Beijing Polytechnic University, 15, 21-26. 

Zhu M. (1990), On the ABS algorithms for perturbed linear systems, Journal of Computer 
Mathematics 29, 181-194. Also as Report 198, Numerical Optimisation Centre, Hatfield 
Polytechnic, 1987. 

Zhu M. (1991a), Finding general solution to perturbed linear systems and its application 
in optimization, Communication on Applied and Computational Mathematics, 5, 61-67. Also 
as QDMSIA 5/90 



37 



Zhu M. (1991b), A note on the condition for convergence of nonlinear ABS method, 
Calcolo 28, 307-314. Also as QDMSIA 17/91. 

Zhu M. and Zhang J. (1995), Solving a least change problem under the weak secant 
condition by the ABS approach, ABSC2. 



38 



